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1 Introduction. 



X 

H Let 



a — (Ai, ... , A„) 
be a list of complex numbers and let 

sk:=\\+ ... +\t fc = 1,2,3,... 

The nonnegative inverse eigenvalue problem (NIEP) asks for necessary and 
sufficient conditions on a in order that it be the spectrum of an entry-wise 
nonnegative matrix. If this occurs, we say that a is realizable, and we call a 
nonnegative matrix A with spectrum a a realizing matrix for a. 

A necessary condition for realizability coming from the Perron-Frobenius 
theorem [2] is that there exists j with Aj real and Aj >| A; |, for all i. Such a 
Xj is called the Perron root of a. 



A more obvious necessary condition is that all the Sk are nonnegative. A 
stronger form of this condition was found independently by Loewy and London 
[11] and Johnson [8], namely: 

(JLL) n " Skm > Sj^, for all positive integers, k and m. 

In terms of n, a complete solution of the NIEP is available only for n < 4. 
The solution for n = 4, expressed in terms of inequalities for the s^, appears 
in the PhD thesis of Meehan [12] and a solution in terms of the coefficients of 
the characteristic polynomial has been published more recently by Torre-Mayo, 
Abril-Raymundo, Alarcia-Estevez, Marijuan, and Pisanero [14]. 

However, the same problem in which we may augment the list a by adding 
an arbitrary number N of zeros was solved by Boyle and Handelman [4] . Using 
a range of tools coming from linear algebra, dynamical systems, ergodic theory, 
and graph theory, they proved the remarkable result that if 

1. a has a Perron clement Ai >| Xj |(all j > 1) and 

2. Sk > for all positive integers k (and Sm = for some m implies s^ = 
for all positive divisors d of m), then 

cr AT := (Ai, ... ,A„,0. ... ,0)(iV zeros) 

is realizable for all sufficiently large N. 

Under these assumptions, a realizing matrix can be chosen to be primitive. See 
Friedland [6] for an extension to the irreducible case. 

The proof of the Boyle-Handelman result is not constructive and does not 
provide a bound on the minimal number A^ — N{a)oi zeros required for realiz- 
ability. 

Finding a constructive proof, with a bound on the minimum numberA'^ of 
zeros required, has been an area of much research, and a number of special cases 
have been resolved. In particular, a best possible result in the case that Re(Aj) 
< 0, for all j > 1, has been obtained by Smigoc and the author [9] and, when 
a is real and has exactly two positive entries, a constructive proof with a bound 
on A^ has also been found [10]. 

In the case that a is real and has just one positive entry, then the inequal- 
ity Si > is necessary and sufficient for realizability. This was proved by 
Suleimanova [13] and this is often viewed as the first result on the NIEP. Fried- 
land [5] re-proved her result by showing that the companion matrix with spec- 
trum a has nonnegative entries, and matrices related to companion matrices are 
used in the cited work with Smigoc. 

Here, a constructive approach to the Boyle-Handelman result is presented. It 
is shown that a certain kind of patterned matrix is " universal" for the realization 
of spectra with power sumss^ > 0, ( fc = 1,2,3, ... ). in the sense that all 
such spectra satisfying the Perron condition (i) above can, with sufficiently many 
zeros added, be realized as the spectrum of a primitive nonnegative matrix with 
that pattern. 



2 A matrix related to Newton's identities 



Let 



Xk :=Mi+ - +M^ fc = 1,2,3, 
q(x) -.^ Ii'^^-^{x - in) 
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The matrix Xn occurs in the context of the Newton identities relating the 
coefhcients of a polynomial to the power sums of its roots. If we use Cramer's 
rule to express the qi in terms of the Xj , we get 

det(X„) = (-l)"n!g„. 

However, the matrix Xn itself, as distinct from its determinant, does not appear 
to have been widely investigated. A key observation is: 



Proposition 1 The characteristic polynomial of Xn is 



Qix) 



- ngix" ^ + n{n - 1)922;" ^ + 



+ n!q„. 



Since Xn has nonnegative entries if the Xi are nonnegative, it follows that the 
spectrum of Q{x) is realizable if the Xi, (i — 1,2, ... ,ri), are nonnegative. 



Suppose that we are given a list a — (Ai, 
as the spectrum of a nonnegative matrix. 
Let 

/(x):=nr.i(x-A,) 
= .t" +pix""^ + 



, A„) that we wish to realize 



Let 

g(.x) :=x" + gix"-i+ ... +g„ 

where 

^^ ^ = 1,2, 



n(n — 1) ... (n — i + 1)' 

Then the corresponding Q{x) is f(x). Now the power sums Xi of the roots 
of q{x) are nonnegative if and only if that holds for 

x"+nqix"^'^+ ... +n"g„. 

Hence we have 

Theorem 2 a is realizable by the matrix X„ if the j th power sum of the roots 
of the polynomial 

2 n—1 

Uf{^)) ■.= X-+PlX--' + ^P2X--' + - ^. r-p3x"-3+ ... + ^L^p, 

n—1 (n— l)(n — 2) (n— 1)! 

is nonnegative for j ~ 1, 2, 3, ... , n. 

But now suppose that we choose N > n and ask for the realizability of a 
with N — n zeros added. This amounts to replacing f{x) by x^^'^f(x) and 
Jnifix)) by x^-"Jw(/(x)), where 

J^(/(x)) := x" + pix"-i + ^J^p2x"-2 (1) 



(A^ - 1)(A^ - 2/^"" •■■ (A^ - 1)(7V - 2)...{N -n+ 1)^"' 

(2) 

So a with N — n zeros added is realizable by the matrix Xn if the j th 
power sums of the roots of the polynomial JAr(/(x)) are nonnegative for j = 
1,2,3, ... ,iV. 

But observe that as A^ — > oo, JN{f{x)) -^ f{x), since n is fixed. 

Suppose that the power sums Sj of the elements of a are positive for all 
J > 1. Then, on continuity grounds, one might expect that for sufficiently large 
N, the power sums of the roots of JN{f{x)) would also be positive. However, 
this is not true in general, but it is true if a has its Perron element 
Ai>|A, I (j = l,2, ... ,n). 

In this case, a with sufficiently many zeros added is the spectrum of a non- 
negative matrix X^. 

Since, we only require that the j th power sum of the roots of J^ifix)) be 
nonnegative for j = 1,2, ... ,iV, one can obtain a bound on the minimal 
number of zeros required. 



Main Theorem 



We now state the main result of this paper. 

Theorem 3 Let 

a — (Ai, ... , A„), 

be a list of complex numbers with corresponding power sums 



Sk := Ai + 



^ni f*^ — 1) ^j "J) ■•■ 



Suppose that 

(i) Ai>| A, I, (allj>l) 

(a) si > 0, and Sm > 0, for all m > 2. 

Let 

/(a;)=.nr=i(x-A,) 



Pn- 



7 = 2max(l,|pi \,\p2 |^/^ ... ,\p„\ 



|1/«N 



and 



Ao =max{| A^ \: j > 1}, 



i? = 



(Ai - Ao) 



3Ai + Aq 
Ai + 3Ao ' 



r = niin(i?, 1), 



max{l, Ai}, A^'o = 



ln(2n - 2) 



N 



\n{i) 
M = min{l,S2, ... ,sno), 



V 3i/2Mr 



Then a with N — n zeros added is the spectrum of the nonnegative matrix 
Xn, with Xk '■= iJ-i + ... + /ij^, fc = 1, 2, 3, ... ,N, where 

Jjv(/(x)) = (,x - N^ii){x - Nfi2)...{X - N^n). 



Given a list a satisfying the hypotheses, it is relatively easy to find N for 
which JN{f{x)) has the corresponding power sums nonnegative, so one obtains 
a reasonably efficient constructive algorithm. However, the number of zeros 
required in the construction is not optimal in general. 



4 Proofs of the results 
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and let C = 
/ 1 

10 

.010 



... 001 

V -n\qn ■ -"^91 / 

be the companion matrix of (5(a;) = a;"+ngia;"^^+n(n — l)(72a;"^^+ •.• + 
n!g„. 

Direct multiplication, using the Newton identities, yields PC = XnP. This 
proves the proposition. 

To obtain the desired bound we use the following refinement by Bhatia, 
Eisner and Krause [3 ] of a classical result of Ostrowski. 

Theorem 4 Let f(x) — x" + aix"^^ + ... + a„ and 
g{x) ^ x"" + bix"-^ + ... +bn 
be real polynomials with roots ai, ... ,a„ and Pi, ... ,/3„, respectively. 
Then there is a labelling of Pi, ... , /?„ such that 



max{|a, -^, \: I < i <n} < i—j=\ (^l^i\ au - b^ 



k I 7 



i—k\llr. 



V3V3/ 
where 7 = 2max{| Cfe \^/^ , \ bk \^'^: l<k<n}. 

[The original Ostrowski result had the factor (2n — 1) in place of ( -^ ]] 



Now let 

f{x) = {x-\i) .... (x-A„) 

and 

g(a;)=a;"+pia;"-i+ ^32;""^+ ... + U„. 

'^^ ' ^ \N-lJ^ \(N-l)...(N-n+l)J^ 

We note that if 5(2;) has nonncgative Newton power sums, then the corre- 
sponding matrix Xj\j is nonnegative and has spectrum A^Ai, ... ,N\n. 

Suppose that Ai >| Xj \ (ah j > 1) and let Aq = max(| \j |: j — 2, ... ,n} 
and R = ^1^. Let £ ^ ^^, so ^ < 1. Let r = min{l,i?}. Let 

Sk = Xi+ - +A^ 
for k = 1,2, ... . Assume that si > 0, and that Sk > 0, for all fc > 1. Let 

Af = min{l,Sfe : k ^ 2,3,...}. 
Let Hi, ... , /i„ be the roots of g{x) ~ and suppose that 
max{| Aj -^j |: j = 1,2, ... ,n}<5, (*) 



where 



with 



Mr 
= 



TO = max{l, Ai}, No 



nNo(m + r)^o-^' 

ln(2(n-l)) 



ln(l/£) 



Then | /xi | is greatest among all the | /Xj |, and , since g{x) has real coefficients, 
Hi is real and, since Ai is positive, so is /ii. Let 

Sk=Hl+ - +^^n■ 

Then \sk-Sk\<E Li I Af - ^f I • Now 

I A,^ - A^- hi A, - M. II A^-i + \t'f^^ + ■■■ + ^^■-' I 

< 6k{Xi + r)''-\ 

Suppose that k > Nq. Then Sk > (Xi - r)'' - (n - l)(Ao + r)'= = (Ai - r)'=(l - 
in - 1)({^)'=) > (^)(Ai - r)'= > 0. For k < No, 

\sk~Sk\<S{l + 2{Xi+r)+ ... +No{Xi+rf"-^) 

< 5nNo{m + r)^""^ = Mr < M. 



So Sk > 0, for all k > 2. Also, Si = si > 0. This shows that if we can 
choose N so that the inequality 

niax{| Xj — /ij |: j = 1,2, ... ,n} < S 

holds for that S, then the corresponding Xj^ will be a nonnegative matrix with 
spectrum Ai, ... , A„ and N — n zeros. Now, 



max{| Xj — ^j |: j = 1, 2, ...,n} < 

v^)(2^fe=i l:Pfc((^_i) ...(iv-fc + i) 



(^)(Ei~.b>-( („_n r^,_t + i) - '>^"" I"" 



But 

TV'^-i , 2n2 

n". 



- 1 < ^^^,if A^ >"2 



(iV-1) ...(7V-fc+l) - N 

By definition, 7 ~ 2max{l, I pfc l^/*^, fc = 1, 2, ... ,n}. Hence 

max{|A,-M, |:j = l,2, ... , n} < (i^)(^)i/". 

But n^/" < 3^/3. Hence 

I67 2^/" 
max{| Aj -Aij |: j = 1,2, ... , n} < — =^^^ < (5, 

provided 

24"+i7" _ 2(167n7Vo(TO + r)^«-i)" 

This gives the required bound. 

There are variations of the Ostrowski bound, some using the Bombieri norm 
in place of the £2 one, available though the work of Beauzamy [1], Galantai and 
Hegedus [7], and these may lead to better bounds for N in certain circum- 
stances. However, the main interest is that such a bound exists, and the general 
form it has. 

When the Perron root Ai = l,a nonnegative matrix A with the given nonzero 
spectrum can be made stochastic. In this case r and £ are measures of the 
spectral gap, which control the rate at which the powers of A converge to the 
stationary state of the corresponding Markov process. The size of A^o is inversely 
related to r and £. 

The number M measures how close to zero the power sums can get, and we 
see its appearance (as M") in the denominator of the bound. 



We conclude with an example involving the realization of a spectrum with 
three nonzero entries. 



Example 5 cr = (/?, cxp(-^),cxp(^j^)) has all its power sums positive if p > 

■y/2 cos(7r/10) = 1.07 // we take p = 1.1, and carry out the algorithm, we 

find that a with 125 zeros added is the spectrum of an 128 x 128 nonnegative 
matrix of the form of X above. The least number of zeros required to be added 
to a to ensure realizability does not appear to be known in this case. 
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